Probing the quantum phase transition in the Dicke model through mechanical vibrations 
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This paper is concerned with quantum dynamics of a system coupled to a critical reservoir. In this context, 
we employ the Dicke model which is known to exhibit a super radiant quantum phase transition (QPT) and 
we allow one of the mirrors to move under a linear restoring force. The electromagnetic field couples to the 
movable mirror though radiation pressure just like in typical optomechanical setups. We show that, in the 
thermodynamical limit, the super-radiant phase induces a classical driving force on the mirror without causing 
decoherence. 
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I. INTRODUCTION 

Quantum phase transitions (QPTs) are qualitative modifi- 
cations that occur in interacting quantum many body systems 
at zero temperature. They are driven by quantum fluctuations 
that are induced by a change either of coupling constants or 
external parameters of the system [lLl2|]- The criticality in the 
phenomena of phase transition is in general governed by a di- 
verging correlation length and, specially in the case of QPT, it 
has been conjectured that a quantum correlation with no clas- 
sical counterpart - the entanglement - must play a crucial role 
Oh- The methods used to study phase transitions in statisti- 
cal mechanics required the introduction of new concepts from 
quantum information theory to characterize and classify the 
critical ground states at the QPT Q]. Seminal articles on spin 
systems undergoing QPT yl S have motivated an enormous 
amount of work aiming at properly defining the role of entan- 
glement in QFTs. For example, the occurrence of some QPT 
are indicated by non-analyticity of entanglement measures 
near the quantum critical point 0,0,0]. However, there is 
no one-to-one correspondence between quantum phase transi- 
tions and the non-analyticity of entanglement measures when 
approaching the quantum critical point Oh- 

Among the most studied systems presenting collective 
quantum phenomena in atomic physics and quantum optics, 
the Dicke model (DM) stands out jlfioll. It comprises a sys- 
tem of 7Y two-level atoms interacting with a single quantized 
electromagnetic mode, and it un derg oes the QFT known as 
super-radiant phase transition Hill — Tl3ll . The finite temperature 
phase transition in the thermodynamic limit of large- N and 
weak coupling of DM was in the seventies fTil [lSl . More 
recently, the role of this QFT on quantum chaos signatures 



has been analyzed 1161 - 2C |, the exact expression for the scaled 
concurrence obtained 12111 . and the scaling of relevant physi- 
cal quantities determined I22I l23h . Universality of these re- 
sults has been extended to the case of inhomogeneous cou- 
pling 03]. On the practical side, some alternative physical 
implementations of effective DM have been proposed in dif- 
ferent settings such as atom-optical systems with multilevel 
atoms and cavity-mediated Raman transitions l25ll . or even 
solid-state devices involving superconducting qubits l26ll . 

Another recent issue related to QPT is the question of how 



a system coupled to a critical bath will behave as the bath 
undergoes a QPT. Topics such as decoherence, the decay of 
Loshmidt echo, and entanglement have been studied in sys- 
tems cou pled to a surrounding environment that presents a 
QPT i27l432ll . In this article, we address this issue of quan- 
tum dynamics of a system coupled to critical reservoir. We 
consider the Dicke model as a structured bath to which a fur- 
ther oscillator is coupled. For this choice, we were inspired 
by typical optomechanical systems where the electromagnetic 
field couples to the mechanical vibrational motion of a mirror 
subjected to a linear restoring force I33l435tl . and we attempt 
to answer the question: What happens to the moving mirror 
when the field and atoms enclosed in the the cavity undergo 
the super radiant phase transition ? We show that this QFT 
induces an effective classical non fluctuating driving force on 
the mirror. This result differs from the common case of a bath 
of harmonic oscillators, for example, where the presence of 
the reservoir causes a. fluctuating driving force responsible for 
the decoherence of the system. 

This article is organized as follows. In Sec. [II] we describe 
the optomechanical system used in this work. In Sec. fXTTJ we 
discuss the thermodynamical limit of this system, and present 
the response of the moving mirror to the super radiant quan- 
tum phase transition. In Sec. IIVI we summarize our results. 
Semiclassical calculations are shown in the Appendix. 



II. MODEL SETUP 

In the Dicke model, an ensemble of N two-level atoms is 
coupled to a single-mode quantized electromagnetic field in a 
cavity. Such a model contains a QFT from the normal phase 
to the super-radiant phase through variation of the strength of 
the coupling constants between the field and the atoms ITTI - 
[I3I1 . The model is described by the Dicke Hamiltonian which 
is given by US (H=l) 



Hr 



where 
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is the angular frequency of the electromagnetic field mode de- 
fined by fixing n (with creation and annihilation operators 
and a, respectively), C the speed of light, L distance between 
the mirrors, loq is the level-splitting of the two-level atoms, 
{J+, J-, J z } the collective atomic operators satisfying the 
usual commutation relations of angular momentum, and 



A = p£ sin(fca;o), 



(3) 



is the coupling constant between the atoms and the cav- 
ity field, with p the atomic electric dipole matrix element, 
£ = sJuj/cqL the r.m.s. vacuum electric field amplitude, eo 
the vacuum electric permittivity constant, and xq the position 
of the atoms in the cavity axis. 

Notice that we have assumed identical coupling constants 
A (same position xq) for all atoms. This means that we are 
considering an atomic gas confined in a region of small extent 
as compared with the cavity size L. This assumption is stan- 
dard in the treatment of the Dicke model lfl9i |2U l27ll . For a 
more general setup where the atoms spread in a considerable 
region of the cavity (gas of large extent), it turns out that the 
critical behavior is essentially the same as that in the ordinary 
model as long as the calculations are done with the average 
coupling constant l24ll . It is worthwhile to mention that the 
generalization for a traveling quantized field presents an even 
richer phase structure when a gas of large extent is considered 
ll36tl . and that the inclusion of atomic motion leads to modified 
normal and superradiant phases ll37ll . 

We now consider that one of the mirrors is allowed to move 
under the action of a linear restoring force, as depicted in 
Fig. Q] This can be achieved, for example, by mounting 
the mirror on a cantilever. The strong coupling regime be- 
tween a mode of the cavity electromagnetic field and the vi- 
brational motion of the mirror has already been experimen- 
tally realized l34ll . Let us denote the angular frequency of 




FIG. 1: (Color online) A collection of N two-level atoms interact 
with a single-mode quantized cavity field. The cavity has one of its 
mirrors fixed and the other is allowed to move under the action of a 
linear restoring force which causes vibration of the moving mirror. 

vibrational motion by io m , and the associated creation and an- 
nihilation operators by c' and c, respectively. This motion 
induces the expected optomechanical radiation-pressure inter- 
action which couples the cavity photon number a' a to the po- 
sition X m = (c' + c)/ y (2mw m ) of the mechanical oscilla- 
tor, where m the effective mirror mass. However, the mirror 
vibration also induces a three-body interaction as previously 
noted for one-atom case HH]. The physical reason behind the 
appearance of this coupling is that once the mirror moves, the 
spatial profile of the cavity field is slightly modified causing 



changes on the strength of the atom-field coupling A. Let us 
now derive such Hamiltonian. Due to the mirror motion, the 
cavity changes it length from L to L + X m . Considering the 
oscillations very small as compared to L, Hamiltonian (Q]) be- 
comes 



nnc ( X m \ f t 

— - II — I a l a + uj m cJc + u} Jz 

p£ 

H — ^[sin(fcxo) - X m 5](a' + o)(J+ + J_), 
v TV 



H = 



(4) 

where k — nir/L and 6 = [sm(kx ) + kxo cos(fcr )]/ L, and 
we also included the free term uj m c i c. We notice that besides 
the radiation-pressure interaction part, a three-body interac- 
tion term also appears. Let us rewrite the above Hamiltonian 

as 



H = 



uja'a + uj m c^c + uJaJ z H — i={<^ + a )(J+ 
v N 



-ga)a{c* + c) j={c^ + c)(a f + a)(J + 

v N 



J-) 

J-), 
(5) 



where g = ujj\L\J (2mw m )] and 77 = p£S/ \J (2muj m ). De- 
pending on the specific choices of the system parameters, 
one or another term may dominate. We choose to work 
with the two-body radiation-pressure term, since it has al- 
ready been experimentally realized fi4\. To compare the rel- 
ative importance of these terms, we need to consider rj/g = 
^[sin(fc:co) + ^ x o cos(fcErj)]. For example, for the funda- 
mental mode n = 1, we can see that rj/g = for kxo ss 
2.02876. This means positioning the atoms around 2L/3. 
With this choice, the system Hamiltonian reduces to 

H = ua) a + oj m c ! c + ujq,J z H — 7=( a ^ + a ){J+ + J-) 

v N 



9° t 1 t ^ 
-— a'a(c' +c), 



(6) 



where we define go = gN. Now the dependence of the ra- 
diation pressure coupling on N, arising from the density of 
atoms p = N / L, appears explicitly in the Hamiltonian. 

III. EFFECTIVE DRIVING FORCE ON THE MIRROR 

We are interested in the thermodynamical limit of this 
Hamiltonian, where the QFT effectively takes place. The 
mathematical procedure to reach this limit in each phase of 
the Dicke model is thoroughly explained in lfl9tl . therefore, 
we shall present here only the points needed to understand 
our model which includes an extra oscillator. 

First, we perform the Holstein-Primakoff transformation 
defined as 



J+ 

Jz 



bWN-tfb 
N 
~2~' 



6 f 6 



(7) 
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where the set of atoms is described by bosonic operators , b, 
satisfying [b, b^} = I. After performing this transformation, 
the system Hamiltonian becomes (up to constant terms) 

H = wa^a + uiob^b + u> m c^c — -j^cfla^ + c) 



Now the system Hamiltonian is in a suitable form for the 
thermodynamic limit to be applied. The Dicke Hamiltonian 
presents two phases, and the physical features of these phases 
motivates the mathematical treatment of the thermodynamical 
limit. In the normal phase A < A c , where A c = y/uJujo/2, the 
simple application of the limit in the Hamiltonian (O leads di- 
rectly to the correct effective Hamiltonian of the normal phase 



H np = H n 



Dickc 



where 



i?° p icke = uja f a + oj tfb + A(a f + a)(& t + b). 



(10) 



(11) 



We can see that in the normal phase, the moving mirror effec- 
tively decouples from the Dicke system and evolves as a free 
harmonic oscillator. We interpret this as a consequence from 
the fact that the Dicke system is only microscopically excited 
in the normal phase. 

On the other hand, for A > A c we must incorporate the 
fact that both the field and the atomic gas acquire macroscopic 
occupations. This motivates the displacement of the bosonic 
modes as al -> a* + yfa and ftt _> fet _ ^ if me 

displacement constants are chosen as 



excited. This is translated as an effective classical driving 
force acting on the mirror. 

This is a very interesting point because an ordinary thermal 
reservoir would never drive coherent oscillations, even at zero 
temperature. There is a special combination of facts leading 
to this particular effect in this optomechanical setup as we are 
now going to explain. It is well known that the coupling to a 
quantum environment results in a fluctuating force acting on 
the system and reflecting the characteristics of the reservoir 
rf39ll . It is precisely the nature of the force fluctuations that 
cause decoherence and damping. On the other hand, if the 
number of elements in the reservoir is small, energy can be 
transfered from the system to the reservoir and then fed back 
to the system (Poincare recurrence) in a period relevant to the 
natural time scales of the problem. In cases where the system 
is coupled to a huge reservoir, the Poincare recurrence time 
will tend to infinity. When thinking in terms of a reservoir- 
induced fluctuating force on the system, the former case cor- 
responds to a practically not-fluctuating force. In the optome- 
chanical system treated here, the mirror is coupled solely to 
a single member of the elements composing the critical reser- 
voir, namely the electromagnetic field. Therefore, Poincare 
recurrences will take place, explaining the energy oscillation 
of the mirror system, i.e., the fact that c^c is no longer a con- 
stant of movement. This still does not explain why the mir- 
ror does not suffer minimal decoherence, something expected 
when the reservoir is traced out, even if it contains only a sin- 
gle element. The point is that the atoms work as a reservoir for 
the field causing it to behave classically, i.e., a — > ^fa such 
that 



50 
N 



a) a(c* + c) — > — aic^ + c) 



N 



(16) 



A 



= -y/N(l-H 2 ) 



y(l-M), 



(12) 
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where /1 = ujujq/AX 2 = (A c /A) 2 , and we set the terms with 
overall powers of N in the denominator to zero after expand- 
ing the square roots in Hamiltonian (0 considering large N, 
we obtain (up to constant terms) 

H srp = H^ + u m c*c-^(l-ii 2 )(c*+cUU) 



where 



H, 



Dickc 
srp 



u> Q {l - /i)(3 + n) 



8/i(l + /i) 



(6 f + bf 



+A/ ViT^ (at + a)(5t + 5) ' (15) 

Therefore, in the super-radiant phase the moving mirror also 
effectively decouples from the Dicke system, but it suffers the 
consequences of the field and atoms being macroscopically 



This is consistent with the fact that the macroscopic excita- 
tion in the field is given by (ahi) /J = a/ J (with 2 J = N) 
in the super radiant phase Il9ll . From the point of view of 
the displacement a* —> + ^/a, it follows that, in the ther- 
modynamical limit, the vacuum fluctuations in the field are 
negligible concerning the interaction with the mirror. 

As an alternative way of showing the macroscopic limit of 
the Dicke system, we present in the appendix A a study of 
the classical analogue of the present model. This method has 
been used to study chaos 11171 [l9|| and the bifurcation of the 
equilibrium in the ground state flal in the Dicke model. It 
gives reliable results in the large N limit of scalable systems 
like the Dicke model. 

From all of this, it follows that the vibrational motion of 
the mirror can then work as a probe to the super radiant phase 
transition. If A < A c , the average occupation (c'c)(t) will be 
time independent according to Hamiltonian ( fTUl ), whereas just 
above the critical point, the mirror starts to evolve under 



with £1 = — :2 ^-(l — /i 2 ), and it will cause the average occu- 
pation of the vibrational levels to change in time. For example, 
if the Dicke system is initially in the ground state of iJ s ^p Ckc 



i.c'c 



+ n( c + + c), 



(17) 
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and the mirror in the vacuum, we will have 



t C )(t) 



2n 2 



[1 - cos(w m t)] 



(18) 



for A > A c . Of course, it is not really necessary to cool down 
the mirror at first place. Any initial state would be good be- 
cause what matters is that c'c is a constant of motion only 
in the normal phase according to our findings. Consequently, 
only in the super radiant phase of the structured bath we will 
be able to see some dynamics of this observable. 

The above behavior in the super-radiant phase is numer- 
ically evidenced in Figf2] for the mean occupation number 
(c^c) (t) simulated with Hamiltonian ([5). The system has been 
initially prepared in the state |^(0)) = | ground) Dicko <8> |0) c , 
where | ground) Dicko and |0) c are the ground states of (Q} and 
H m = ui m c)c, respectively. Although far from the thermody- 
namic limit, we can already see that as J increases, the oscil- 
lations tend to the regular energy oscillation ( fl~8l of the forced 
harmonic oscillator. The simulation with finite J was per- 
formed using standard diagonalization packages. When map- 
ping the atomic Pauli operators to collective angular momen- 
tum operators, it is easy to see that the basis for treating the 
atomic part in the Dicke system has actually dimension 2 J+l, 
where J = N/2. Consequently, cases with J about 15 are 
amenable to implementation in current personal computers. 
There is here, of course, the complication that the basis for 
the whole system includes also the field and the mirror sub- 
systems. However, since we are simulating just the cases in 
which the initial state comprises the fundamental state of the 
Dicke system and the vacuum state of the mirror, the numer- 
ical calculations converge quickly without the need to drasti- 
cally increase the basis dimension. 




under the presence of the Dicke system (reservoir). In Figf3] 
we keep the same initial preparation considered in Figf2j and 
we show the numerically evaluated time evolution of the Von 
Neumann entropy S c {t) = — Tr[p c (t) lnp c (t)] of the mirror 
state p c {t) = Tr a ,/p(t), where p{t) is the density operator 
of the whole system. One can clearly see the tendency of the 
entropy to diminish as J increases, confirming the statement 
above (S C (Q) = for initial state of the mirror). 




FIG. 3: (Color online) Time evolution of the Von Neumann entropy 
of the mirror state for increasing values of J. The system parameters 
are again uj — ujq = 10io m and A = 0.6 (super radiant phase). 



IV. CONCLUSIONS AND FINAL REMARKS 

We have presented a new effect in optomechanics, which is 
a super-radiance generated classical driving force on a moving 
mirror. We started from the Dicke Hamiltonian and derived an 
effective Hamiltonian considering harmonic motion of one of 
the mirrors composing the cavity. The subsequent application 
of the thermodynamical limit reveals the appearance of a driv- 
ing force on the moving mirror in the super-radiant phase of 
the Dicke system, without causing any coherence loss during 
the time evolution of the mirror state. Although the thermody- 
namical limit of this system may not be completely under the 
reach of present technology, we expect that the tendency of 
the mirror to appear driven in one phase and free in the other 
may soon be observed in the regime of few atoms in optome- 
chanical systems of the type reported in l40ll . 



FIG. 2: (Color online) Time evolution of mean occupation number 
(c*c)(t) for increasing values of J. It is clear the the curves tend to 
the driven case d 1 8t found in the thermodynamical limit (TL). The 
system parameters are u = ujq = 10oj m and A = 0.6 (super radiant 
phase). 

According to our findings, the reservoir does not induce de- 
coherence on the mirror in the thermodynamical limit. This 
can be seem from the fact that there are no products of a or 
b with c in Hamiltonians ( fTOb and ( TBI . Consequently, the 
initial amount of entropy of the mirror state will not change 
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Appendix A: A classical analogue of the model 

In this appendix, we present a classical analogue of the 
model using a time dependent variational method based on co- 
herent states ll4lll . This method produces for the normalized 
coherent states representation the following classical Hamil- 
tonian (in units of K = 1 = m): 



H = (z a ,w,z c \H\z a ,w,z c ) 

= y(9i +Pi) - + ^(q 2 + P2) + — (9 3 



2 A 



9l92 



,9o 



V2 



4.7 



Pi) 

(Al) 



where is given by ©, Hi = (qf + p?)/2, and the co- 
herent states |js a , to, z c ) are defined as |z a ,m, z c ) = \z a ) <E) 

\w) ® |z c ) where = e ( 2dt ^|0)(d = a,c) and |w) = 
\J,—J). The phase space coordinates (<?i,pi) 



(1+ww) 

and (q2,p 2 ) (similarly for (53,^3) with z c ) read |4 



9i 
Pi 
92 

P2 



1 + Ml 



(til + 



J (w — w) 
1 + ww i 

(z a + Z Q ), 

1 fca ~ Z~a) 

2 i ' 



(A2) 



(A3) 



According to (lAU . the associated classical equations of 
motion are given by 



A 



9i 
Pi 
92 



w pi - 
-w gi 
ujp 2 



Pi 9i 92 



2J\/2J-H! 
2A 



y/2J-H 



A 



2J 
2A 



P293 



p 2 = -^<72 - 

i>3 = -u m q-$ 



y/2J-n! 



iq-2 



91 



9?92 



V2J V 2 J - Hi 
.9oV2_ 



2.7 



"9293 



.9o%/2 2 2 
^j-(92 + 



(A4) 



In the macroscopic limit (N = 2 J — >• 00), and below crit- 
ical coupling A < A c = where no bifurcation of the 
minimum energy fixed point occurs, one can see from the last 
two equations that the oscillator with frequency oj m decouples 
from the rest of the system and follows free evolution. 



On the other hand, at A = A c a pitchfork bifurcation oc- 
curs for the Dicke model lowest energy fixed point lfl8ll . such 
that the stable fixed points for A > A c are (<7in,Pio) = 



(TV2^(l-^),0) and (q 20 , P20) = (±a/4J^(1 - /x 2 ), 0). 
As a kind of first order perturbation theory, we may use this 
values into the last two equations in (1A4I) . and this leads nat- 
urally to a forced harmonic oscillator equation in the limit 

N -> 00: 



93 + ^m93 



0, ifA < A r 



V2u m go£(l- V 2 ), ifA>A c 



(A5) 



Now, it is clear that, at least for small go-coupling, the driving 
term has the same form as indicated in Sec. [Till i.e. propor- 
tional toil = — ga£p(l — /J 2 ) and consistent with the quantum 
Hamiltonian ( fT71 ). 



The Dicke model quantum phase transition considering dis- 
sipation in the a cavity field has been treated in ref. l25ll . 
where it has been shown that the transition point is shifted 

to A c = ^\Jljuio(1 + ^|), where k is the loss rate. More- 
over, the fixed point of semiclassical equations results in 

71 [920 



W2Q\ 



= ±V2J^-Jl- u 2 . For small 

v OJ—ZK V t 

losses (k -C u>), the right hand side of above Eq. jA5i can be 

substituted by s/2uj m go^(l — M 2 )(l — t)- Since the mean 
photon number in the super radiant ground state is written as 

z a Za = r?20 " 2 ~ P20 and this is diminished for k > 0, the force 
acting on the mirror in the super radiant phase is also dimin- 
ished due to dissipation. 
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